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The Analogy between Lesage’s Theory of Gravitation and the 
Repulsion of Light. 

By G. H. Darwin, F.R.S., Plumian Professor and Fellow of Trinity College, 
in the University of Cambridge. 

(Received May 2,—Read May 18, 1905.) 

I am not aware that anyone has taken* the trouble to work out Lesage’s 
theory, except in the case where the particles of gross matter, subjected to 
the bombardment of ultramundane corpuscles, are at a distance apart which 
is a large multiple of the linear dimensions of either of them. Some years 
ago I had the curiosity to investigate the case where the particles are near 
together, and having been reminded of my work by reading Professor 
Poynting’s paper on the pressure of radiation,* I have thought it might be 
worth while to publish my solution, together with some recent additions 
thereto. 

If a corpuscle of mass m moving with velocity v impinges on a plane 
surface, so that the inclination of its direction of motion before impact to the 
normal to the surface is $, it communicates to the surface normal momentum 
kmv cos S, and tangential momentum k'mv sin &; where k is 1 for complete 
inelasticity, and 2 for perfect elasticity, and k' is 0 for perfect smoothness and 
1 for perfect roughness. 

In the following paper the effects are investigated of the bombardment by 
Lesagian corpuscles of two spheres, which are taken to be types of the atoms 
or molecules of gross matter. The effects of the normal and tangential 
components of the momentum communicated by each blow from a corpuscle 
will be treated separately. 

§ 1. The Normal Component of the Impacts. 

Suppose that there are in space n corpuscles per unit volume, moving 
indiscriminately in all directions with velocity v; let the mass of each 
corpuscle be m, and p the density of the medium, so that p = mn. 

Suppose, further, that the corpuscles are so small that collisions between 
them are rare enough to be negligible. The system thus described is that of 
Lesage’s mechanism for explaining gravitation, provided that the corpuscles 
are not perfectly elastic. 

* ‘ Phil. Trans.,’ A, vol. 202, pp. 525—552. I have had the advantage of corresponding 
with him on the present subject. 
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388 Prof. G. B. Darwin. Analogy between Lesage’s [May 2, 

As stated in the introductory remarks, when any corpuscle strikes a plane 
surface at such an angle that S is the inclination of its motion before impact 
to the normal to the plane, it communicates to the plane normal momentum 
kmv cos $, where h is a factor lying between 1 and 2, being 1 for perfectly 
inelastic corpuscles, and 2 for perfectly elastic ones. 

If two bodies, A and B, are exposed to the bombardment of corpuscles, the 
Lesagian attraction between them will vanish if the corpuscles are perfectly 
elastic, because the corpuscles reflected between A and B will on the whole 
cause such a repulsion between the bodies as exactly to counterbalance the 
attraction. 

But in the present investigation we shall only consider those corpuscles 
which come from infinity, and pay no attention to their behaviour after 
impact on. the surfaces of A and B. This can, I suppose, be only completely 
justified when the corpuscles are perfectly inelastic. 

If the velocity of each of the n corpuscles in unit volume be represented 
by a vector drawn from a centre, the ends of the vectors will lie equally 
spaced on the surface of a sphere of radius v. If from the centre of the 
sphere of velocities a cone of small solid angle Soo be drawn, the number 
of corpuscles per unit volume, whose directions of motion fall within the cone 
is nSco/4:7r. 

Let us consider the pressure on a surface exposed to bombardment. On the 
surface draw a unit area, which may clearly be treated as plane, and let $ be 
the zenith distance of the axis of the cone 8co. On the unit area draw an 
oblique prism, with edges of length v, parallel to the axis of the cone 8co; the 
volume of this prism is v cos $. Hence the number of corpuscles, whose 
directions of motion fall within the cone 8co, lying inside the oblique prism 
at any moment is nv cos S 8co/4:7r. At the end of unit time elapsing after the 
moment under consideration all these corpuscles will have struck the plane, 
and each will have communicated momentum towards the nadir equal to 
kmv cos 3. Thus the momentum towards the nadir communicated by this 
class of corpuscles is kmnv 2 cos 2 $Sa)/47r. 

Since mn = p, the pressure (or momentum towards the nadir communicated 
per unit time and area) is 

-- £ v — 2 cos 2 $So), 

47 T 

where the summation is effected for all directions whence corpuscles may 
come. 

If the particles may come from all over the sky, the summation is for half of 
angular space, and2 cos 2 S8co — §7r. Therefore the pressure^ on the surface 

is § iven lj y P = \i P v\ 
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1905.] Theory of Gravitation and the Repulsion of Light. 389 

If we put k = 2, this is the well-known result of the kinetic theory of gases. 
If, however, a certain cone of space is screened off, so that no corpuscle can 
come from thence, the pressure is given by 

p = ±kpv 2 —-ff - S cos 2 SSco, 

4:77“ 

where the summation is taken over the screening a,rea. 

In fig. 1 let the origin be at the unit area, and let P be the intersection with 
the celestial sphere of the normal to the area; let the screen be circular 

Fig. 1. 


z 



round the Z axis, and let AA' be one quadrant of the screen. Let the radius 
of the screen be a, and the zenith distance (for an observer standing on the 
unit area) of the centre of the screen be 6. Then 

ZA = ZA' = ZP = ft. 

Let 0 , <£ be the polar co-ordinates of an element sin 0 d0 d<fi of the screen 
at Q, so that the Sco of the formula for the pressure is sin 0 d0 d<f>. Also, since 
the 5 of the formula is the angle PQ, we have 

cos 5 = cos 0 cos ft ~f~ sin 0 sin /3 cos <j>. 

Hence, taking the summation over the area of the screen, 

1 1 r 27r / a 

~ £ cos 2 3Sco = — j J (cos 0 cos ft + sin 0 sin ft cos <fi) 2 sin 0 d0 dcj> 

= j* (cos 2 0 cos 2 ft +1 sin 2 0 sin 2 ft) sin 0 d0 
Jo 

= [ [(f cos 2 ft —J) cos 2 0 -f- \ sin 2 /3] sin 0 d0 

J o 

= \ (cos 2 ft cos a sin 2 a + cos 3 a — cos a). (1) 

Thus 


p = ±kpv 2 — ±Jepv 2 (cos 2 ft cos asin 2 a + § -f-^cos 3 a—cos a). 

2 e 2 
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390 Prof. G. H. Darwin. Analogy between Lesages [May 2, 

]STow let there be two spheres exposed to bombardment, a sphere A with 
radius a , and a sphere B with radius b ; and let B be the distance between their 
centres. As stated in the introduction, A and B are meant to typify atoms or 
molecules of gross matter. 

Let an element at P of the sphere A be determined by colatitude 9 from the 
axis BA, and by longitude <fi measured from some fixed plane passing through 
BA, and let the element be expressed by a 2 sin 0 d6 defy. 

This element is screened by the sphere B , and therefore the screen is 
circular. The semi-angle of the screening cone was denoted above by a, and 
the observer’s zenith distance of the axis of the screening cone was /3. 

Fig. 2. 


x 



If then we denote the distance PB from the element to the centre of the 
sphere B by r , we have from fig. 2, 

b - n B, - a o Bcos6 — a 

sm « — - , sin p = — sin 6, cos p =-, 

t r r 

r 2 — B 2 + a 2 —2Ba cos 6. 

It is to be noted that 6 has not the same meaning as that used in 
evaluating (1). 

If the angle QPA is obtuse, the whole of the sphere B is visible to the 
observer at P. This case will be considered first, and afterwards those regions 
on the sphere A from which B is partially below the horizon of the observer. 

We want to evaluate the total resultant force on the sphere A, due to the 
pressure on all its elements, in the direction A to B. As far as concerns the 
element under consideration, the resultant force in the required direction is 
—p cos 0 . a 2 sin 0 d9 d<f>. In order to determine the whole of that portion of 
the resultant which is due to pressure on the portions of A from which B is 
completely visible, this expression will have to be integrated from 0 = 2tt to 
<f) = 0, and from cos 0 = (a + &)/B to cos 0 = 1—that is to say, from the 
colatitude in which B just begins to set until B is in the zenith. 

The constant part of the expression for the pressure, namely, \hpv 2 , is the 
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1905.] Theory of Gravitation and the Repulsion of Light. 391 


same all round the sphere, and it cannot contribute anything to the resultant 
force; hence it may be omitted. 

Since the expression for the pressure does not involve </>, the integration 
with respect to </> simply introduces the factor 2i r. Thus the resultant force 
estimated towards B, in as far as it is due to the elements from which B is 
wholly visible, is 


j[cos 2 /3 cos a sin 2 a, -f1 d -g- cos 3 a —cos a] cos 6 sin 6 d6 , 


where the limits of 6 are cos 1 ( 


and 0. 


This expression may be integrated rigorously, but it would not be expedient 
to do so, because the result is very complicated, and because I have failed to 
obtain a rigorous integral for the remaining portion of the complete expression, 
namely, where B is only partially visible. 

It will now save trouble if we omit the factor §7 rkpaPv 2 and reintroduce it 
later, and if besides we choose B as our unit of length. 

Hence we have 


b . n 1 . /) Q cos 6 —a 
sm a = -, sin B — - sm 6 , cos B — - 

r r r 


1 + a 2 — 2a cos 6, 


and the limits of 6 are cos” 1 (a + b) to 0. 

Cos a can be expanded convergently in powers of 5/r, so long as the two 
spheres do not actually touch, and we easily obtain 

f-^COS 3 «-COSa = 1—+ - 1 L_ + 5 3 r _+..... 


also 


sin 2 ol cos a 


_ } W 

ry >2 2 ^,4 8 1 IT ^,8 


cos 2 /3 


• (cos 2 6 — 2a cos 6 -f a 2 ). 


Hence the integral becomes 

\— (cos 2 0 — 2a cos 6 + a 2 ) cos 6 sin 6 d6 

+ W\(^i + %~ r 5+t\~s---) cos 6 sin 6dO. 

Now r 2 = (1 + a 2 ) [ l — ■§ cos ■ 


If, then, we write 

A 1 


we have 


2a B b 2n 

~ 1 + a 2 ’ 2n (l + a 2 ) n+1> 

r 2 = (14- a 2 ) (1 — A\ cos 6 ). 
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Whence 



Integrating between the appropriate limits, we find 



} 






It will be observed that the suffixes of the i?’s and of A\ indicate their 
orders in the powers of 1 JB. 

This formula may be applied to each term of the preceding integrals, so 
that the result may be obtained in terms of A x and B 2n . 

The several terms must then be expanded in powers of a 2 , and terms 
rearranged in their several orders. 

After some tedious analysis I find, on reintroducing B, which was 
treated as being the unit of length, the following result as far as the ninth 
order:— 



1 + JL JL W) H a(8a' L + Ma 2 b 2 + 2lb*) 


3.5 B 3.5.7 B* 3.5.7A 5 




WV^ 3 ) . 
( 2 ) 



This formula gives the force on the sphere A towards B in so far as it is 
produced by the normal component of the bombardment on those parts of A 
from which B is completely visible. It remains to investigate the regions 
whence B is only partially visible. 

In fig. 3, as in fig. 1, let the origin be at a unit area on sphere A, and let P 
be the intersection of the normal to the area with the celestial sphere. Then 
if PC is a right-angle, CY will be the horizon which bounds the screen. Let 
Z be the axis of the screening cone, and DBE the cone itself. As in the 
previous case, the semi-angle of the cone is denoted by a, and PZ is /?. 

The half of the screen which is shown in the figure is CBEZ. 

If Q be any point in the screen, the pressure on the unit area on the 


sphere A is 


,2 


,2 
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Fig. 3. 



the summation being carried out over the whole screen. As before, we may 
neglect the constant part of the expression for the pressure. 

In this figure more than half the sphere B is visible, and the case when less 
than half is visible will be treated subsequently. 

If 0 , <j) are the polar co-ordinates of Q, 

-!- 2 cos 2 PQ 8co = 2 ”||( cos ^ cos ft + sin 0 sin ft cos <f>) 2 sin 0 d0 d<f>, 

where the integral is taken over the area CBEZ and the other half not 
shown in the figure. 

The area of integration may conveniently be divided into two parts—viz., 
the part inside the dotted circle and the part between CF and BE. 

The part of the result corresponding to the dotted circle may be obtained 
at once from our previous result (1), for CZ the radius of the cone is equal to 
\ir—ft, and therefore the result is 

{ cos 4 /3 sin /3 +1 + i sin 3 ft —sin ft }. 


For the remaining portion, take the limits of integration so as only to 
include CBEF; then doubling the preceding formula, so as to permit us to 
halve the range of </>, we write the integral in the form 



0 cos 2 ft + \ sin 2 0 sin 2 ft + \ sin 2 0 sin 2 /3 cos 2 <f> 

+ 2 sin 0 cos 0 sin /3 cos /3 cos <£) sin 0 d0 d<f>. 


If B is any point on the boundary CB in colatitude 0 , its longitude cf> is 
equal to 7 r—CZB ; also CZ = \n r—/3. Then from the right-angled spherical 
triangle CZB we have 

cos cf) = — cot /3 cot 0. 


Thus the limits of </> are cos 1 [ — cot ft cot 0] to 0 ; and the limits of 0 are 
a to \it — ft. 
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394 Prof. G. H. Darwin. Analogy between Lesage’s [May 2, 


At the upper limit of (f> we have 

sin (f> = ,/(l-cot 2 /3cot 2 0) = 


and 


sin /3 sin 0 

*0.1 COS iS COS 8 // ' o a o /o\ 

sin 2cf> = -2 -tL- v/(sin 2 0-cos 2 /3). 

sin J /3 sm J 6/ 


Integrating with respect to <£, we obtain 

-y|(cos 2 0 cos 2 /3 + | sin 2 0 sin 2 /3) cos -1 (—cot /3cot (9) d8 

— —f cos /3 cos 0 ^(sin 2 8 — cos 2 /3) sin 8 dO 
27tj 

+ -1cos /3 cos 0 y/'( sin 2 8 — cos 2 /?) sin 6 d6 . 

The last two terms involve the same integral and only differ in the coeffi¬ 
cients ; they therefore fuse together with a coefficient 3/27r. 

The whole may now be written 

— - (cos -1 (— cot /3cot 8)[} cos 3 8 (J- cos 2 /3— -|) -f1- sin 2 /3 cos 0] dd 
7 rJ do 

+ — (sin 0 cos 0 cos /3 ^/(sin 2 0 — cos 2 /3) c?0. 

27rJ 

The last term is a perfect integral, and the first may be integrated by 
parts. Effecting this, I find that the indefinite integral of the whole is 

__ JL cos 8 [cos 2 8 (cos 2 /3 — -^)-h sin 2 /3] cos -1 (— cot /3 cot 8 ) -f 

+ G -^£ (sin 2 8 - i) ^(sin 2 8 - cos 2 j3). 

Zi7T 

Then taking the expression, between limits, and adding the part corre¬ 
sponding to the dotted circle, we have 

\ sin ft [sin 2 /3 (cos 2 /?—i) + sin 2 /?] 

—— cos a [cos 2 cl (cos 2 /3—£) + sin 2 /3] cos -1 ( — cot (3 cot a) 

2ir 

,1 . cos 0 1 

4-— sin 1 -r—- — - 
o 7 r sin u o 

+ C ° S -^ (sin 2 a —-*-) y/ (sin 2 a — cos 2 /3) 

2nr 


+ \ [cos 4 /3 sin /3 + f 4 - - 5 - sin 3 —sin /3]. 

This expression admits of reduction, and it gives 

— S cos 2 PQ So)= i—i cos# [f + sin 2 a(^--cos 2 /3)]cos -1 (— cot/3 cot 8) 

2ir 0 27t 

■f ~ cos /9 (sin 2 a — j)\/(sin 2 « — cos 2 /3)-+ -5- sin -1< ??--. (3) 
2ir o 7 r sm a 
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The pressure is the same quantity multiplied by J Tcpv 2 
In the case when less than half of the sphere B is visible, the dotted circle 
vanishes. It appears that at the upper limit 6 = a at the lower 6 = /3-~|-7t; 
the angle cos" 1 ( — cot /3 cot 6) is the angle (f> at the point C, and therefore 


vanishes ; and at the same point sin 1 is equal to — \ir. It is, then 


easily verified that the formula (3) holds good, and we may therefore apply it 
both when more and when less than half of B is visible. 

If p is the pressure as given by this formula, and a? sin 6 d6 d(f> is the 
element of area of the sphere A, the component of pressure in the direction 
AB is —jp cos 0. a 2 sin 0 d# d</>. This has to be integrated over the 
“ sunset ” region of the sphere A. The integration with respect to is 
effected at once, and we are left with 



cos 1 (a + b), when B is taken as the unit of length. 

We have accordingly to multiply (3) by 7r/rpaVsin 6 cos 6 dd and integrate 
it from 6 = cos" 1 (a — b) to cos" 1 (a-b&). 

I have not been able to effect the integration rigorously, and shall not give 
the full details of the tedious analysis involved in the approximate integration. 
It is convenient, in the first place, to change the independent variable. 

Putting then sin \[r = ? ~ , so that the limits become + \i r, we have 




If we drop the factor nrhpa?v 2 temporarily, it will be found that the subject 
of integration may be written in the form 


i& 4 (1 -fi 2a?) (£ + sin 2 yjr) — b 6 (sin 2 ^— f) + - cdfi yfr sin y/r -f- sin 2 ijr) 

7T 1 

+ J— ab h sin 2 \[r cos yfr (2 sin 2 ylr + 13). 

07T 

— a6 5 sin (J-f-sin 2 -x/r)—1 + 2<x 2 ) (^-f sin 2 ^)—yjp ( s in 2 

Air 


— a6 5 sin (J-f-sin 2 yfr)——^fr{b 4: (l + 2a?) (^-f sin 2 ^)—(sin 2 yjr — ^)} 


— ty{l + 2a 2 ) sin yfr cos (2 sin 2 ^-fi 13) 
b 6 

+ — sin yfr cos yfr (£ + § sin 2 yfr — sin 4 tIt), 

07T 


multiplied by (ab cos ty-V 2 sin ^ cos yfr) d\jr. 
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The first two lines of the first factor are even functions of y/r, and the last 
three lines are odd functions; the first term of the second factor is even, and 
the second term is odd. Then since the limits of ^ are r, it follows that 
we need only multiply the first two lines of the first factor by ab cos tyd-fr 
and the last three lines by — 6 2 sin cos ^ All the integrals involved 
are then known, viz.:— 


. 2 
sin 2w i Jr cos yfr dyfr = 

I-frr 


2^+1 9 


r&r . 1 

yfrsmyjr COS yfr d\jr = - 7 r, 

J -*«■ 4 

f iw g 

yfr sin 3 COS yjr dyfr == —7T, 

-i*r 


32 


f" Sill 2 " 1 yjr COS 2 ™ ~ tjffi~ ^ ^ • 1 " ■> ^ 

J_ ijr r r r (2m + 2n)(2m + 2n— 2)... 2 

By means of these I find that the component force on the sphere A from 
A towards JB due to this pressure is 

^ a+ +^a« 3 + ¥^+^ 2 -tW& 3 )}- 


Now on comparing this with the expression for the first part of the force 
in (2) we see that it exactly annuls the terms of the sixth and eighth orders; 
and we may feel confident that the * term of the tenth order would be 
similarly annulled. The result is of so simple a character that it must 
surely be possible to prove it in some shorter way. However this may be, 
our final result for the resultant force on A towards B , which I may call 
¥ a , is given by 


?A = 



8a, 4a(4a 2 -f 7b 2 ) a(8<x 4 + 36& 2 6 2 -f 21& 4 ) 

15B + 105 ^ + 105 ^ 


+^M^ 6 +^ 2 +^+W)}- (4) 

We may observe that, as a rough approximation, 


?A = 


i-rrkp 


v 2 a 2 b 2 

(k-tVO 2 ’ 


By symmetry the full expression for the force on B towards A is 



4& (46 2 -f 7& 2 ) 
105^ + 


b(8¥ + 36a 2 b 2 +2la i ) 
105B b + 



In considering the excess of F^ above Fg, it will suffice if we drop the 
terms of the seventh and ninth orders, since the reader will easily be able to 
extend the result if he desires to do so. 
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Then F a -F b = frrkp . * 



]■ 


It follows, therefore, that, as the result of the normal component of the 
bombardment, the larger sphere is more strongly urged towards the smaller 
than the smaller towards the larger. 

| 2. The Tangential Component of the Impacts; Effect of loth Components 


together . 


When a corpuscle m strikes a plane surface so that $ is the inclination of 
its motion before impact to the normal to the plane, the tangential component 
of its momentum is mv sin $, and it communicates to the plane tangential 


momentum equal to h'mv sin $, where h' lies between 0 and 1, being 1 for 


perfect roughness and 0 for perfect smoothness of the surface. 

Following a procedure exactly similar to that adopted previously, we see 
that the corpuscles whose direction of motion before impact lie in a small 
cone of solid angle communicate a tangential force to the plane equal to 


sin $ cos $ Bo). 


47T 


Let Z, m, n be the direction cosines of the normal to the plane bombarded, 
and V, m', n' the direction cosines of the axis of the cone So). Then the 
direction cosines of the projection of the axis of the cone on the plane, which 
is identical with but opposite in direction to that of the tangential force 
communicated to the plane, are 


V — ZcoSrS m'—mcos$ 
sin $ ’ sin 5 


n' — n cos S 
sin S 


Hence if X , Y, Z are the components of, the force thus imparted to the 
plane by the bombardment, 


X — — — h‘ r pv 2 X cos $ {V — l cos -9) 8oo, Y = — j- h'pv^cosS^m'— mcos$) S®, 

47r 47T 

Z — — -i- k'pv 2 H cos $ (n' — n cos Sco, 

47T 


where the summations are carried out over all the directions from whence 
corpuscles arrive at the plane. 

If corpuscles come from all over the sky, X — Y == Z = 0; hence the 
forces are 



where the summations are carried out over the screening area. 
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In fig. 1 the direction cosines of the normal to the bombarded surface are 
l = sin ft, m — 0, n— cos ft ; 


and V — sin 8 cos <f>, m' — sin 8 sin <£, n = cos 8 ; 


cos 3 = cos 8 cos ft + sin 8 sin ft cos <£. 


The limits of integration are 


$ =27t to 0, 8 = ol to 0, and Sco = sin 8 dd d<f>. 


Therefore 


2 r2ir ra 

X = — k'pv 2 [sin 6 cos <£ (cos 6 cos ft + sin 8 sin ft cos cf>) 

^ r7r Jo Jo 

— sin ft (cos 6 cos ft + sin 8 sin ft cos<£) 2 ] sin 8 d8 dcf>, 




Integrating with respect to </>, we find 



= — \k'pv 2 sin ft cos 2 ft sin 2 « cos «, 


Y = 0, 

Z = Ik'pv 2 sin 2 ft cos ft f (f cos 2 0—-J) sin 8 dd 


o 


= \k'pv 2 sin 2 ft cos ft sin 2 ol cos a. 


In fig. 2 P is at the bombarded element of surface, and is the origin of 
fig. 1; then we take PB as the axis of Z , and the dotted line PX as the axis 
of X. Accordingly the component force in the direction A to B is 

Z cos (ft — 8) + X sin (ft— 8) = \k f pv 2 sin ft cos ft sin 2 ol cos ol sin 8. (5) 

The element of surface at P is a 2 sin 8 d8 d<f>, and it is to be noted that 8 is 
used in a different sense from that employed in evaluating X, Y, Z. 

Hence the total force on the sphere A from A to B, due to the tangential 
component arising from the bombardment, for all those parts of A from which 
B is wholly visible, is 


\k'pv 2 a 2 sin ft cos ft sin 2 a cos a sin 8 . sin 8 dd dcf>, 


the limits of <£ being 2tt to 0, and of 8 from cos -1 (a + b)jR to 0. 

The integration for <£ involves only multiplication by 2ir. 

As before, I take temporarily B as unit of length, so that the limits of 
8 are cos _1 (a + b) to 0. 








Downloaded from rspa.royalsocietypublishing.org 


o 
z 
oc 


2<Z}U 


^u 
Sq.o3i/> 



Ss 


uj U 

EO 

H 


C/5 

O 



u 

o 

0 * 

CL 


_r a 
< z 

y 5 

b m 

2 <ziu 

SS=!v5^ 

h^zS 
5>; uj — 
Si v u 

Sq.o 3 to 



£ £ 


w(j 

KO 

H 



1905.] Theory of Gravitation and the Repulsion of Light. 399 


Then 


b 2 , 5 4 


sm a cos a = - 2 -§^-t^..... 
Hence the integral becomes 


sin ft cos ft = sin 0 (cos 0 — a). 


|tt /;>- 2 « 2 | (JJ - i _ i . j sin® 6 (cos 6~a)d0. 

If we use the same abbreviations as before in the development, in inverse 
powers of r, and note that 


we find 
[ b 2 ‘ 




sin 3 9 cos m+1 6dQ = 


(i a -f b) m+2 (a -f &) m+4 

--T“- 


(m+2)(m + 4) to + 2 m + 4 


f b 2n -i^^co" — Bi f ^ ^ ^ (%+l)(w + 2 ) 2 , 

]^+ 2 Sm t ' cosy ^-^\ 2 . 4 + !, O 1 + - 2 l-476 1 


(w+l)(w + 2)(?i + 3) 2 . 3 

+--- __. Al + ... 

<a + bf (« + l) (a + bf A (n+l)(n + 2) (a + 5 ) 4 A 2 
2 lT 3“~ 1 Ti-4~ 1 


I ■O + fr ) 4 ■■ («- + ! ) (a+ff A (w + l)(w + 2 ) («+5 ) 6 ; 2 
41! 5 2! 6 

f ^ 2 ” sin 8 ^# — Bi ^ ^ ^ 41 . ( n +l)( n + 2) 2 4 0 , 

J r 2» + 2 sm ^^-^2»| 0 +- rr 0 ^i+-2l-375 1 + - 


(a + 6 ) _(m,+ 1 ) (a + b) 2 A (n+l)(n + 2) (a + bf A 

1 1 ! 2 1 2 ! — 3 


(a+5) 8 (w + 1 ) («+5) 4 . (» + !)(»+2) (a+ 5)6 . 2 

3 1! 4 1-r 2 ! 5 ‘ 1 

These formulae may be applied to each term of the integral, and the result 
is thus obtained in terms of ^4i and B 2n . We next expand these in powers of 
a 2 , and rearrange the result in their several orders. Finally on reintroducing 
R I find the following result:— 

1 7 / v 2 a 2 b 2 8 a 4<x(4a+7&) &( 8& 4 + 36& 2 & 2 +21& 4 ) 

^ 9 R 2 \ 3.5 R 3.5.7 R 3 3.5.7.fi 5 

- J 7 (7W aS + 2 8 t« 4 ^ 4- f + e 6 ®) 

~^g>2-j>i(% a2 +* ah - V 62 )-^ (|a 4 + -V«a 3 5 + x y ( 2 52 _ s| a j3 + j,_ 

(6) 
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400 Prof. G. H. Darwin. Analogy between Lesage’s [May 2, 

This gives the force on the sphere A towards B as arising from the 
tangential component of the impacts, in so far as it is due to those parts of 
A from which B is completely visible. It remains to consider the regions 
whence B is only partially visible. 

By proceeding in the same way as before I find, for the portion of the 
result which does not involve that part of the screen which is represented by 
the dotted circle in fig. 3, after the integration with respect to </>, 

X = - k'pv 2 [ [sin/3 cos 2 /3 (£— f cos 2 0) cos -1 ( — cot/3 cot 0) 

2tF V-i3 

-f cot /3 (-J — f sin 2 6) cos ^^/(sin 2 0— cos 2 /3)] sin 6 cl6. 

Z = k'pv 2 [ [sin 2 /3 cos /3 (f cos 2 6 — J) cos -1 ( — cot 6 cot 6) 

2lT J^TT-ZS 

4 (1 — f cos 2 /3) cos 0^/(sin 2 6 — cos 2 /3)] sin 6 d6. 

The last terms in each of these expressions are integrable as they stand, 
and for the first terms of each 

j*(| cos 2 6 —J) cos '" 1 ( — cot /3 cot 6) d0 

= \ sin 2 6 cos 6 cos ” 1 (cot /3 cot 6) 4 -cos /S^/^sin 2 6 — cos 2 /3). 

In proceeding to the limits, it is to be noted that cos"^-—cot/3 cot 6) is n r, 
when 0 = -| 7 r—/3. Thus after integration we shall have, for that portion 
which depends on the term involving cos -1 ( — cot /3 cot 0), when 0 = \ir — /3, 

X = 4 k'pv 2 .sin 2 /3 cos 4 /3, Z = — k'pv 2 . ^ sin 3 /3 cos 3 /3. 

Air A Air A 

How, since the resultant force is given by ^cos (/3—0)4 JTsin (/3 — 0 ), this 

portion of the resultant is equal to — ^ k'pv 2 sin 2 /3 cos 3 /3 sin 0, where 0 has 

the meaning indicated in fig. 3. 

The force due to the portion of the screen represented by the dotted circle 
is given by the result (5), namely \k'pv 2 sin /3 cos /3 sin 2 a cos a sin 0, when a 
is put equal to ZC, which is equal to | 7 r—/ 3 ; it is, therefore, equal to 
Ik'pv 2 sin 2 /3 cos 3 /3 sin 0, and exactly annuls the term referred to above as 
resulting from the term in cos -1 (—cot /3 cot 0) with 0 = \ir— /3. Thus I 
find that the complete values for X and Z ’ inclusive of the dotted circle, are 
given by 

X = ~ k'pv 2 { — sin /3 cos 2 /3 sin 2 a cos a cos " 1 (—cot /3 cot a) 

4iir 

+ -i- cot /3 [sin 2 a—cos 2 /3 — 3 sin 2 a sin 2 /3]y/(sin 2 a—cos 2 /3)} 

Z = -i- k'pv 2 {4 sin 2 /3 cos /3 sin 2 a cos a cos -1 ( — cot /3 cot a) 

Xrr 

—•-i- [sin 2 a — cos 2 /3 —3 sin 2 a sin 2 /3]v/(sin 2 a — cos 2 /3)}. 
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These formulae may be shown to be equally true of the case where less 
than half of B is visible, and they are therefore applicable throughout. The 
component force of A towards B is equal to Z eos (/3— 6) + AT sin (/3 — 9), and 
becomes 



In this formula 6 is used in the sense indicated by fig. 3. 

This expression has to be multiplied by a 2 sin 6 d6 d(f), and integrated from 
<j> = 27t to 0, and from 6 = cos -1 {a — b)jR to cos -1 (a + b)/R. The integration 


with respect to </> merely involves multiplication by 2i r. ."For the integration 


with respect to 6, I change the variable to yjr, and, as before, develop the 
expression and effect the various integrations. 

The final outcome of some tedious analysis is that the result is the same as 
the last three terms of ( 6 ) with the sign changed. Hence, when we add this 
contribution to the force to ( 6 ), those three terms simply disappear. 

T^hus the tangential component of the impacts give as a resultant, say 
¥'a, acting on the sphere A towards B } 


y' — i i' v 2 a 2 b 2 f h 8 a 4&'(4&+76) a ( 8& 4 + 36a 2 6 2 + 21b 4 ) 

A “ s' 9 JR* l 3.5 5 ■ 3.5.7 & 3XW 



a*b 2 + ±a%* + ±¥) k ( 7 ) 


As a rough approximation we have 


If we form F'j, or the force acting on the other sphere, it is clear that 
F 'a *— F ' B has a form similar to that found previously for F^ — F 5 , but 
it has the opposite sign. Hence, as the result of the tangential component of 
the bombardment, the larger sphere is less strongly urged towards the 
smaller one, than the smaller towards the larger. 

On comparing (7) with the result for F^ in (4), we see that they only 
differ in the signs of all the terms after the first, and in the fact that k' 
replaces k . The result is of so simple a character that it is probable that it 
may be derived by some elementary considerations which escape me. 

Thus, including both the tangential and normal components, we have 


F^ + F^. 





(8> 
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402 Prof. G. H. Darwin. Analogy between Lesages [May 2, 

In the case where the momentum of the impinging corpuscles is completely 

absorbed by the surface struck, we have k = k' = 1 , and the force is 

v 2 a %2 
^7Tp—j^ . 

Thus, in this case the force varies rigorously as the inverse square of the 
distance. No doubt this simple result may be proved much more shortly in 
the case of the complete absorption of the momentum of impacts.* In any 
other case the result can only be regarded as approximate, because we have 
neglected reflected particles. But in general it seems certain that the inter¬ 
action between the two spheres will not be equal and opposite. 

§ 3. Repulsion of Radiation. 

When the Lesagian corpuscles deliver their whole momentum on impact, 
we have the exact converse of the case of radiation, for in the one case we 
consider all the particles which converge on to a given element of surface, 
and in the other case they all diverge. A corpuscular theory of light would 
give the same result as the electro-magnetic theory as regards repulsion ; hence 
we see that two radiating and perfectly absorbing spheres at the same 
temperature will repel one another rigorously as the inverse square of the 
distance.! The case of a perfectly reflecting sphere which receives radiation 

# Such a proof is given in a note to the following section. 

t [June 23, 1905.] Professor Larmor has given me a direct proof of the above result. 
This I paraphrase in my own words as follows :— 

A sphere radiating from its surface in the manner of a perfect radiator may be replaced 
by a uniform distribution of radiating spherules inside it, if we suppose the radiations 
from the several spherules not to interfere with one another. This follows from the fact 
that the radiation issuing towards any zenith distance is proportional to the depth of the 
crowd of spherules beneath it in the given direction. Now for a sphere that depth is a 
chord of the sphere, and is therefore proportional to the cosine of the zenith distance. 
Hence the law of radiation of the crowd of spherules is the same as the natural law of 
radiation from the surface itself. 

Suppose that one of the spherules P emits n' corpuscles of mass m with velocity v per 
unit time, and let it be distant r from the centre A of an absorbing sphere of radius a. 
With origin at the spherule measure colatitude 6 from PA, and longitude from some fixed 
plane passing through PA. 

The number of particles emitted per unit time through solid angle sin 6 dd d<p is 
n' sin 6 dS dcfrl^n. Each of them carries momentum mv, and therefore the component of 
momentum along PA absorbed per unit time by the sphere A is 

n'mv 

sind cos 6d6d<t>. 

To find the whole thrust on the sphere A we must integrate thus from <f) — 2rr to 0, 
and from 6 — sin -1 ~ to 0. 
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1905.] Theory of Gravitation and the Repulsion of Light. 403 

is analogous to the case treated in § 1 , where we consider only the effect of 
the normal component of the impacts. If two perfectly absorbing spheres 
have different temperatures, the action and reaction between them will not be 
equal and opposite. The analogue of this case in the Lesagian hypothesis 
would be that the velocities of the particles which strike one of the spheres* 
should be different from those which strike the other. 

If one side of a body of any shape be at a higher temperature than the 
other side, the body will be subject to a force tending to propel the cooler 
side forward, and to drag the warmer side after it. This follows from the fact 
that the recoil of the emission on the warmer side is greater than that on the 
cooler side. The result that a hot sphere will pursue a cold one is a special 
instance of this more general conclusion. 

In the Lesagian hypothesis the pressure has a definite relationship to the 
amount of energy received per unit area and per unit time, and the law: is 
the same as that which governs the relationship between radiation and the 
recoil of light. The result has so great a physical interest that it seems 


Hence the resultant repulsion exercised by the spherule on the sphere is 

o? 

irimv - 2 • 

Now suppose that there are N spherules arranged uniformly in a sphere B , of radius £>, 
whose centre is distant R from the centre of the sphere A. 

Since each spherule repels the sphere A inversely as the square of the distance from 
its centre, it follows, as in the theory of attractions, that the aggregate of them repels 
inversely as the square of the distance between the centres of the spheres A and B. 

Hence the total repulsion of the crowd of spherules must be equal to 

i , 7 , a? 

It remains to find the value of Nn\ in terms of the equivalent radiation from the 
surface of the sphere B. 

Nn' is the total number of corpuscles emitted per unit time by the whole crowd, and 
this must be equal to the total number emitted from the surface of the sphere B. 

Now we have taken in the text above n to represent the number of corpuscles emitted 
from unit area of the surface of the sphere. Hence we have 

Nn r = 47 rb 2 n. 


It follows that the repulsion between the spheres is 

o?b 2 

Trnmv • 

But it appears from the latter portion of § 3 that 

nm = \<jv = %pv. 

tfa 2 }) 2 

Thus the repulsion is J np ? and this is identical with the result obtained at the 


end of § 2 above. 

VOL. LXXVI.—A. 
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worth while to investigate this matter more closety. We have already seen 
in § 1 that the number of Lesagian corpuscles, which strike unit area in 


unit time is % — ^° - S - 8co , and each corpuscle before impact carries energy 


\m r v 2 . Hence if the surface completely absorbs all the energy on impact, we 
can at once find /, the total absorption of energy. It is given by 



But nm is the mass of the Lesagian medium per unit volume, and has been 
denoted by p, and therefore 


Complete absorption of energy corresponds to the case h = 1, and there¬ 
fore the Lesagian pressure is p = j?pv 2 . 



Hence 


This result is the converse of the case of a corpuscular theory of radiation, 
and I is then the radiation of the surface, whilst p is the pressure of radiation. 

We shall now see how this same result may be obtained, when the 
subject is considered from the point of view of radiation. 

Suppose that n cos 98ooj r iT be the number of corpuscles of mass m emitted 
per unit time with velocity v from unit area of a surface through an 
elementary cone Sco in zenith distance 9. The total number of corpuscles 
emitted per unit time is n, because 



Each corpuscle carries ^mv 2 energy, and the radiation is given by 

I = \mnv 2 . 


But mn has not the same meaning that it had in § 1, and therefore we must 
consider what it is. 

The number of particles radiated per unit area and time through the 


Tb COS 9 

elementary cone Sco towards zenith distance 9 is - Sco, and therefore 


7T 


for a small element of area Ss it is n CQS ^ Sco 8s. 


IT 


We may then concentrate the radiation from 8s at its centre, and consider 
the distribution of corpuscles emitted. Since the corpuscles have mass m 
and move with velocity v, the total mass of the particles at any moment in 


an element of volume r 2 8co distant r from the centre of 8s is 


cos 98co8s. 
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Hence the density of corpuscles radiated from 8s at a point distant r from it, 
in zenith distance 6, is mn CQ ^ -- 8s. 

7TVT J 

We may now find the density at a point distant z from an infinite plane 
radiating surface. Take origin in the plane vertically under the point at 
which the density is to be found, and let p, <j> be the distance and azimuth of 
an element of radiating surface. Tor such an element 


8s = pSpSd), cos 6 — - . 77 —% ——, r 2 = p 2 + z 2 . 

\/(p 2 +z 2 ) 


Then the density at the point under consideration is the sum of the 
contributions of all elements of the plane, and if we denote that density 
by <7, we have 



Hence 


The same result will be true infinitely near a curved radiating surface. 

We found above in considering Lesagian bombardment that I = ^pv s . 
Hence it appears that cr = ±p. This result might have been foreseen, because 
near the plane half the corpuscles are screened off, and so the density must 
be half that in free space. 

The normal recoil of the radiation is clearly equal to the normal!component 
of the momentum radiated per unit time and area, and therefore 



Completing the integration and substituting for mn its value, we have 


The density of energy in space is clearly \<rv 2 , and this is equal to \p 
or 2 Ijv. 

In the electro-magnetic theory of light the density of energy in space is 
.also 2 Ijv, but I understand that the pressure is only half that computed from 
the corpuscular theory and is equal to f Ijv ; thus in that theory the density 
■of energy in space is three times the pressure. 

§ 4. The Resistance to a Sphere moving with Uniform Velocity. 

We may find the force on such a sphere by imparting to the Lesagian 
corpuscles a uniform drift in superposition on their common velocities v. 

Let us first find the pressure on a plane unit area. Take the normal to the 
urea as axis of and suppose that the direction cosines of the axis of an 


2 f 2 
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406 Prof. G. H. Darwin. Analogy between Lesages [May 2 , 

infinitesimal cone Sco are sin 0 cos <fr , sin 0 sin <£, cos 0. Then, if U, V ’ W, are 
the components of the uniform drift of the corpuscles, the velocity of the 
corpuscles whose direction of motion is parallel to the axis of the cone is 

v + U sin 0 cos cf> + V sin 0 sin cf> + W cos 0. 

The momentum towards the nadir communicated by this class of corpuscles 
per unit time is therefore 

(y + U sin 0 cos cf) + V sin 0 sin <£+ W cos cj>) 2 cos 2 0 Sco. 

We must then take Sco = sin 0 d0 dcfi>, and integrate through half of angular 
space. The integration with respect to </> may be effected at once, and£the 
expression becomes 

\kp [*\v 2 + i ( U 2 + V 2 ) sin 2 0+W 2 cos 2 0 + 2vW cos 0] cos 2 0 d0 

J o 

= \kp (V + K ^ 2 + V 2 )(l— p?) + W 2 p 2 + 2vWfi\p 2 d f j J 
Jo 

= ^p[v 2 + i(U 2 + V 2 + W 2 )+fW 2 +fvW]. 

Now, suppose a sphere of radius a to be moving with velocity u in the 
direction of the axis of Z, and that a 2 sin 0 d0 dcj> is any element of its 
surface. 

When we reduce the sphere to rest, we have the case just considered, and 
W of the preceding investigation is equal to u cos 0. The portion of the 
pressure corresponding to the terms v 2 + \{L 2 + V 2 + W 2 ) is the same 
all round the sphere, and will produce no effect; it may, therefore, be 
omitted. 

Then the component of the pressure on the element in the direction of the 
new axis of Z is 

jrkp [| u 2 cos 2 0 + J- mi cos 0] cos 0 . ct 2 sin 0 d0 defy. 

This must be integrated all over the sphere. After integration with 
respect to <f>, which merely involves multiplication by 2ir , we have 

\irkpa 2 f (f u 2 cos 2 0 + vu cos 0) cos 0 sin 0 d0 

J o 

= \irkpo? j* vup?dfji — -j7rJtpct 2 vu. 

If the whole momentum is absorbed, Jc = 1. 

Next consider the tangential component of the momentum. In the investi¬ 
gation of § 2 , put 1=0, m = 0 , n = 1 ; V = sin 0 cos cf), m' = sin 0 sin <£, 
n' = cos 0. Then the components of force are 
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X = — 4 k'phs in 0 cos 0 cos $ \v 4- U sin 6cos<f>+V sin 0 sin <£+ W cos 

47r 

Y = . . . . . . sin <£ . . .. . 

0 . 

These must be integrated through half of angular space. Effecting the 
integration with respect to <j>, we have 

X = —\k'p f^sin 9 cos 9 (v U sin 9 + 77JF sin 0 cos 6) sin 6 d6 , 

Jo 

F = —\h'p j^sin 6 cos 6 (v Fsin 6 4* VW sin 9 cos 9) sin 9 d9. 

J o 

From this we easily find 

X = - tk'p (v U+ fe UW), Y = - ik'p (v F+ T 8 ^ VW). 

In applying this to the case of the sphere, we have, as before, W = u cos 9 , 
and if we take the meridian as the X axis of the preceding investigation 
U = — u sin 9, and V = 0. Hence 

X = £k'p (vu sin 0 + feu 2 sin 0 cos 0). 

The component force on the sphere, parallel to the motion u, is — X sin 9. 
Hence, the whole force on the sphere, due to the tangential component, is 

pa 2 ||(w sin 9 -f -feu 2 sin 9 cos 9) sin 2 9 d9 dcj> 

j -^TT 

vu sin 3 9 d9 = \vrk! pd 2 vu. 
o 

For complete absorption of momentum h' = 1. Hence, in this case the 
normal and tangential components contribute an equal amount to the 
resistance, which becomes 

§ irpa 2 vu . 

The uniform normal pressure, say p, on the sphere when at rest is %pv 2 . 
Hence, the resistance is 

4:7 ra 2 p - . 
v 

I take it that a similar formula gives the resistance to a radiating sphere, 
because the recoil corresponding to the emission of a particle is exactly equal 
to the momentum communicated by one impinging. 

§ 5. Summary. 

Various hypotheses may be adopted as to the form and constitution of the 
elementary portions of matter and of Lesage’s ultramundane corpuscles. 
In this paper I consider the elementary portion of matter to be a sphere, 
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408 Prof. Gr. H. Darwin. Analogy between Lesages [May 2, 

and I suppose the sphere to be either smooth or rough, and the corpuscles 
to be either perfectly elastic or inelastic in their collisions. 

With perfectly smooth spheres and perfectly elastic corpuscles, it is clear 
that the total energy of the system remains unchanged.when two spheres, 
immersed in the medium, are made to approach or recede from one another. 
As no work is done by such movements, there can be no force on either 
sphere. It would be excessively difficult to prove the vanishing of the force 
between the two spheres by a detailed examination of the impacts - of 
the corpuscles, because it would be necessary to take into account the 
corpuscles which are reflected from either sphere so as to strike the other. 
It is, however, certain that in this case there would be no force, and there¬ 
fore imperfect elasticity or roughness in the spheres are necessary conditions* 
for the applicability of Lesage’s theory. 

In the case of partial elasticity and roughness, it would be even more 
difficult than in the former case to trace the effect of reflected corpuscles 
which strike the other sphere. But in proportion as the inelasticity and 
roughness increase, so will a solution, which only takes into account first 
impacts, increase in accuracy. We may fairly conjecture that a very moderate 
degree of inelasticity would suffice to make such a solution fairly correct. 
However this may be, no attempt is made in this paper to consider these 
repeated impacts. 

A fundamental objection to the physical truth of Lesage’s hypothesis 
lies in the fact that it demands a continual creation of energy at infinity to 
supply the gravific machinery. But Lord Kelvin has suggested a manner 
in which this physical absurdity might be avoided.* He supposes that the 
corpuscles are capable of absorbing energy in the form of internal agitation. 
On each impact some of the energy of translation is converted into energy 
of agitation, and a repartition of the energies of translation and agitation is 
effected by the mutual collisions of corpuscles according to Clausius’s law. 
If, however, the work of this paper is correct, this suggestion will not serve 
to remove all the defects of Lesage’s hypothesis. 

I here suppose two spheres to be subjected to bombardment, and evaluate 
the effects of the normal and tangential components of the several impacts 
separately. It is thus possible to make various hypotheses as to the degrees 
of elasticity and roughness. It appears that neither the normal nor the 
tangential components of the impacts give rise to forces of attraction 
between the spheres which vary rigorously as the inverse square of the 
distance between their centres. In fact the resultant force acting on one 

* “On the Ultramundane Corpuscules of Lesage,” £ Phil. Mag.,’ May, 1873, vol. 45, 
p. 321, fourth series. 
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of the spheres, due to the normal component, varies (approximately) as 
the inverse square of a distance equal to the distance between the two 
centres diminished byof the radius of the sphere in question; and the 
resultant force, due to the tangential component, similarly demands the 
augmentation of the distance by of the radius of the sphere. 

If the two spheres are unequal in size, then, as far as concerns the normal 
component, the diminution of the distance, so as to maintain the law of 
inverse square, is greater for the larger sphere than for the smaller; and 
the like is true for the augmentation in the case of the effect of the tangential 
component. 

It follows that if the normal component is only effective, the larger sphere 
is more strongly impelled towards the smaller than the smaller to the larger, 
and the converse is true for the separate action of the tangential component. 
In general, the sum of the two effects will not insure equality of action and 
reaction, nor the rigorous truth of the law of the inverse square. If these he 
necessary conditions for the truth of any theory of gravitation, then Lesage’s 
hypothesis and Lord Kelvin’s modification stand condemned. It is true that 
the inequality of action and interaction may be avoided by supposing that 
all elementary portions of matter are rigorously of the same size, but this still 
leaves the law of inverse square imperfectly fulfilled. 

There is, however, one limiting case in which these particular imperfections 
in the theory are avoided. If the inelasticity of the corpuscles is complete 
and the roughness of the sphere such as absolutely to annul the tangential 
velocity of a corpuscle on impact—in other words, if the absorption of 
momentum on impact is total—the law of inverse square becomes rigorous, 
and action and reaction become equal. This supposition leaves the necessity 
for the creation of energy at infinity in its acutest form. 

The case of the total absorption of energy on impact is strictly analogous 
to the repulsion of light, for the emission of light may be regarded as the 
exact converse of Lesage’s mechanism. Thus the preceding investigation 
proves that two radiating and completely absorbing spheres at the same 
temperature repel one another rigorously as the inverse square of the 
distance between their centres. 

If they are not at the same temperature they will tend to move (as indi¬ 
cated by Poynting) so that the cooler sphere leads. This appears to be a 
special case of a more general law, namely, that a body with one end hot and 
the other cold will tend to move with the cold end leading, because the recoil 
of the emission of radiation from the hot end will be greater than that from 
the cold end. 

Another effect of the recoil of the emission of light may, perhaps, be of 
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410 Lesages Theory of Gravitation and the Repulsion of Light . 

importance in solar physics. Poynting shows that the impulse of solar 
radiation at the earth is 5*8 x 10“ 5 dyne per square centimetre. Now the 
earth is distant 214*4 solar radii from the sun; hence the recoil of light at 
the sun's surface must be 5*8 x 10“ 5 x (214*4^ 2 = 2f dynes per square centi¬ 
metre. An “ atmosphere ” is about 10 6 dynes per square centimetre, and 
although the pressure just computed is very small compared with a terrestrial 
standard atmosphere, yet its effect may be worthy of consideration. 

In the last section the resistance is evaluated which an isolated sphere 
suffers when moving with uniform velocity in a Lesagian gravific medium. 
It appears that the resistance is equal to the area of surface of the sphere 
multiplied by the pressure per unit area on a surface at rest, and by the ratio 
of the velocity of the sphere to the velocity of the gravific corpuscles. I take 
it that the same result will give the resistance to motion of a radiating 
sphere. 
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